function [c,ceq] = fminconstr(x,eC,ell,theta,sigma,pi,rho,eta,i,M,A,B)

c = []; % no nonlinear inequality
ceq = sys(x,eC,ell,theta,sigma,pi,rho,eta,i,M,A,B); % the fsolve objective is fmincon constraints